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This paper considers the role of acoustic impedance, flow resistance, configuration and geometrical 
dimensions of regenerators on the power produced in travelling-wave thermoacoustic engines. The 
effects are modelled assuming a pure travelling-wave and ideal gas, which allows defining a pair of 
dimensionless factors based on the “net” acoustic power production. Based on the analysis provided, 
the acoustic power flow in the regenerators is investigated numerically. It is shown that impedance 
essentially reflects the proportion between the acoustic power produced from heat energy through the 
thermoacoustic processes and the acoustic power dissipated by viscous and thermal-relaxation effects 
in the regenerators. Viscous resistance of the regenerator mainly determines the magnitude of the volu¬ 
metric velocity and then affects the magnitude of acoustic impedance. High impedance and high volu¬ 
metric velocity are both required in the regenerators for high power engines. The results also show 
that the optimum transverse dimension of the gas passage exists, but depends on the local acoustic 
impedance. In principle, it is possible to obtain an optimum combination between these two parameters. 

© 2009 Elsevier Ltd. All rights reserved. 


1. Introduction 

In broad terms, thermoacoustic effect refers to the energy trans¬ 
fer between a compressible fluid and a solid material in the pres¬ 
ence of an acoustic wave. There are typically two practical 
implementations of the thermoacoustic effect: imposing an acous¬ 
tic wave in the compressible fluid along a solid material leads to 
hydrodynamic heat pumping effects (due to the compression and 
displacement of the fluid elements), which in turn result in gener¬ 
ation of a temperature gradient within the solid material in the 
direction of the sound wave propagation. Conversely, imposing 
an appreciable temperature gradient within the solid material 
may result in spontaneous generation of an acoustic wave along 
the direction of temperature gradient in its vicinity. These two 
interactions form the basis for engineering thermoacoustic coolers 
(or heat pumps) and thermoacoustic engines, respectively. 

The thermoacoustic effect was first qualitatively explained by 
Rayleigh, who formulated the following criterion [1 ]: “If heat be gi¬ 
ven to the air at the moment of greatest condensation, or be taken 
from it at the moment of greatest rarefaction, the vibration is 
encouraged”. However, it was not until 1960s when a much more 
accurate quantitative theory was developed by Rott [2,3], who de¬ 
rived the wave equation and energy equation for a single frequency 
sound propagating along a temperature gradient in a channel. 
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There are two basic designs of thermoacoustic engines imple¬ 
mented so far: standing-wave and travelling-wave. In the stand¬ 
ing-wave engines [4-6], a thermoacoustic core (a thermoacoustic 
“stack” sandwiched between the hot and cold heat exchangers) is 
placed in a half-wavelength or quarter-wavelength standing-wave 
resonator. Imperfect heat transfer between the gas and the solid 
surface is used to introduce a significant time delay between the 
gas motion and the gas thermal expansion (or contraction) to meet 
Rayleigh’s criterion. By analogy with RLC electrical oscillators that 
dissipate much electrical power to sustain their oscillations, the 
standing-wave engines are fundamentally inefficient due to the 
imperfect heat transfer. The highest efficiency, defined as the ratio 
of acoustic power running out from the ambient heat exchanger 
over the heat input into the hot heat exchanger, reported in the 
open literature is only 18% [6]. Limited by the intrinsically irrevers¬ 
ible thermodynamic cycle, the thermal efficiencies of standing- 
wave engines are thought to have an upper bound of about 20% [7]. 

From the point of view of thermodynamics, to improve the ther¬ 
mal efficiency of standing-wave engines, one has to improve the 
heat transfer between the gas and the solid surface. However, this 
contradicts the basic design feature of the standing-wave engines: 
that is their reliance on the imperfect heat transfer to achieve the 
required time delay to sustain the oscillation. Ceperley [8,9] was 
first to realize that when a travelling sound wave propagates 
through the regenerator (where the thermal contact between the 
gas and solid material is excellent) from the cold to the hot end, 
the heat transfer interaction between the gas and the solid 
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material undergoes a Stirling-like thermodynamic cycle: the gas 
experiences thermal expansion during the displacement toward 
the higher temperature and thermal contraction during the dis¬ 
placement toward the lower temperature (see Fig. la). In this 
way, the correct time phasing is obtained to meet the Rayleigh’s 


criterion to produce power. Strictly speaking, however, the ther¬ 
moacoustic core (a regenerator sandwiched by the hot and cold 
heat exchangers) just works as a power amplifier. To complete this 
Stirling-like thermodynamic process, acoustic power has to be fed 
to the ambient end of the regenerator with a near travelling-wave 
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Fig. 1. Schematic diagrams of different types of travelling-wave thermoacoustic engines. The arrows show the directions of the acoustic power flow; (a) schematic of the 
thermodynamic cycle that a “gas parcel” experiences in the regenerator; (b) one-wavelength looped tube travelling-wave thermoacoustic engine; (c) bypass-type travelling- 
wave thermoacoustic engines; (d) torus-type travelling-wave thermoacoustic engine; (e) hybrid geometry with a travelling feedback. 
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phasing. Therefore, either an external acoustic source, or a feed¬ 
back mechanism should be introduced. For example, Sugita et al. 
[10] used an electro-dynamic acoustic source to test the travel¬ 
ling-wave power amplifier for cryo-cooler applications. In Gardner 
and Swift’s [11] cascade thermoacoustic engine, one stage of a 
standing-wave thermoacoustic core is used to produce acoustic 
power, which is fed into the adjacent travelling-wave stages for 
power amplification. Various acoustic feedback mechanisms were 
used in travelling-wave engine designs by Yazaki et al. [12], de 
Blok [13], and Backhaus and Swift [7,14]. 

Yazaki et al. [12] were probably the first to demonstrate a prac¬ 
tical travelling-wave thermoacoustic engine, which is schemati¬ 
cally shown as Fig. lb. Here, the thermoacoustic core was placed 
in a one-wavelength “looped tube”, which transports part of the 
acoustic power (out from the hot end of the thermoacoustic core) 
back to the ambient end. However, this engine proved very low 
efficiency, as later realized due to the very low acoustic impedance 
of the working gas in the regenerator, which caused large viscous 
losses resulting from high acoustic velocities in the regenerator 
and the resonator [12]. At almost the same time, two other differ¬ 
ent travelling-wave engines were proposed independently: de Blok 
[13] demonstrated a travelling-wave engine with a bypass feed¬ 
back mechanism, which is schematically shown as Fig. lc. Back¬ 
haus and Swift [7,14] proposed a different type, in which 
thermoacoustic core is placed within a torus with a length much 
shorter than the acoustic wavelength at the operating frequency. 
A long standing-wave resonator is connected to this torus just after 
the secondary ambient heat exchanger. It is schematically shown 
as Fig. Id. The thermal efficiency of the latter (defined as the acous¬ 
tic power delivered to the long resonator over the heat input) is so 
far the highest reported: 30%, equivalent to 41% of the theoretical 
Carnot efficiency. 

Essentially, the torus (without the long resonator) of Backhaus 
and Swift’s travelling-wave engine is an acoustic analogue of the 
free-piston Stirling engine. The main idea contained in their inno¬ 
vation is to mimic the free-piston Stirling engine using a compact 
acoustic network (i.e. the torus): using acoustic inertance, compli¬ 
ance and resistance. Its high efficiency was achieved by several 
new means: firstly, they used an “inert” feedback (inertance tube) 
to bring the timing phase close to zero in the regenerator. By anal¬ 
ogy with RLC electrical oscillators, this mechanism added a posi¬ 
tive feedback to the thermoacoustic power amplifier, thus, the 
self-oscillation can be sustained. Secondly, the compact network 
obtained a high acoustic impedance in the regenerator to suppress 
the dissipation: namely 15-30 times of p m a/A, where p m is the 
mean density of the gas, a is sound speed and A is the cross-sec¬ 
tional area. Thirdly, the thermal buffer tube (TBT) and the “jet 
pump” were utilized to suppress the Rayleigh streaming and Gede- 
on streaming which caused the parasitical heat losses. Finally, a 
long quarter-wavelength branch resonator was introduced to the 
torus, which provides the acoustic resonance to let the torus work 
at a relatively low frequency (compared with one-wavelength 
mode of the torus). As a result, the lengths of the parts in the torus 
are much smaller than the wavelength, which reduced the acoustic 
dissipation in the torus. 

However, it should be noted that, the role of this long stand¬ 
ing-wave resonator is as harmful as it is helpful. In essence, it is 
a Helmholtz acoustic resonator, its large surface dissipates much 
of the acoustic power delivered from the torus. If this dissipation 
is counted, the achievable thermal efficiency (useful acoustic 
power over the heat input) will be less than 0.30. In essence, 
this standing-wave resonator has two functions at the same 
time: a resonator to achieve the low frequency resonance for 
the torus and a parasitic (harmful) acoustic load. Theoretically, 
one should attempt to design a resonator which can sustain 
the resonance but with low harmful dissipation. Several 


researchers tried different shapes and sizes of resonators to 
achieve this goal [15-17]. It was reported that the resonators 
of alternative form can improve the engine’s performance in 
terms of the pressure amplitude. However, no improved thermal 
efficiency has been reported so far. 

An alternative to using acoustic resonators is to introduce 
mechanical resonators into the thermoacoustic engine design. 
Sugita et al. [10] presented such an approach which used a pair 
of solid displacers as the resonator instead of a long resonator tube 
to realize both the standing- and travelling-wave thermoacoustic 
pressure generators. A self-sustained acoustic oscillation was ob¬ 
tained in their engine. Here, the displacers are functionally similar 
to those in free-piston Stirling engines. However, this solution is 
not practical due to the high cost of low-dissipation mass-spring 
oscillators (displacers mounted on flexure bearing), and due to 
adding moving parts, which were to be avoided in thermoacoustic 
engines in the first place. 

An alternative strategy to reducing the harmful parasitic dissi¬ 
pation is to design a mechanism which can provide the resonance, 
but at the same time directly extract the acoustic power. Backhaus 
et al. [18] demonstrated such a travelling-wave thermoacoustic 
electric generator. Instead of the long resonator, they connected a 
pair of moving coil linear alternators to the torus. The mechanical 
resonance of the linear alternators induced and sustained an 
acoustic resonance in the torus, while the alternators utilized the 
acoustic power directly to produce electricity. After a careful opti¬ 
misation, the engine efficiency (acoustic power to the alternator 
pistons over the heat input) of 24% and thermal-to-electric conver¬ 
sion efficiency (electricity power to heat input) of 18% were 
obtained. 

The most recent development on the travelling-wave engines is 
the work by de Blok [19]. He pointed out that in the torus-type 
travelling-wave engine (see Fig. Id), the combination of the high 
regenerator impedance and the large acoustic loss in the stand¬ 
ing-wave resonator make the onset temperature difference very 
high. To reduce the onset temperature difference and apply this 
technology to utilize solar and waste heat, he proposed hybrid con¬ 
figuration with a travelling-wave feedback waveguide, which is 
schematically shown as Fig. le. 

In fact, this new configuration can be looked upon as an im¬ 
proved variation of Yazaki’s looped engine: (i) the phase shifters 
(inertance tube and compliance) are introduced to tailor the 
acoustic phasing to a near travelling-wave phasing at the cold 
end of the regenerator; (ii) a bypass is added to shunt part of 
the volumetric velocity, that is to turn it away from the engine 
core branch. In this way, a travelling-wave looped resonator 
(consisting of the feedback wave guide and the bypass tube) is 
achieved. In essence, the acoustic wave keeps running in this 
looped resonator with a near travelling-wave phase almost 
everywhere. The engine core branch extracts part of acoustic 
power from this looped resonator, and then tailors the acoustic 
phase to a near travelling-wave phasing and feeds it into cold 
end of the regenerator. After the amplification in the regenera¬ 
tor, the amplified acoustic power is injected back to the looped 
resonator. Because of the near travelling-wave phase, the acous¬ 
tic loss in the looped resonator is quite low. Furthermore, the 
regenerator was deliberately designed with relatively low 
impedance and a large cross-sectional area. Both of these mea¬ 
sures reduced the dissipation in the regenerator. As a result, 
the onset temperature difference was reduced significantly. The 
average onset temperature difference was reduced to 65 K, 
which is the lowest reported so far in the open literature. 

The other advantage of de Blok’s engine is that the “soft spot” 
(near travelling-wave phasing region) for the regenerator is rela¬ 
tively long, thus, multiple (i.e. multistage) thermoacoustic cores 
can be utilized in one engine. This would be very helpful for 
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increasing the power density when using low-temperature heat 
sources. However, although the relatively low regenerator 
impedance helps to obtain the low onset temperature difference, 
it is necessary to consider the regenerator impedance when apply¬ 
ing this type of engine in high pressure amplitude situations when 
a high power density is required. On the other hand, the very short 
regenerator is a particular design for utilizing the low temperature 
difference. If one wants to apply this design to high temperature 
heat source to get high power density, the heat loss from the hot 
heat exchangers to the ambient heat exchangers (especially in 
multistage design) will be considerable. Unfortunately, so far, there 
is insufficient experimental data to prove that this new type of en¬ 
gine has a better efficiency than the torus type with a standing- 
wave resonator. 

In parallel to the work on devising an appropriate acoustic net¬ 
work (the resonator designs discussed above), there have also been 
a number of studies focused on the optimisation of the regenera¬ 
tors for travelling-wave engines. Backhaus and Swift [20] fabri¬ 
cated and tested a parallel-plate regenerator, which has the 
dimensions similar to the tested stainless stacked-screen regener¬ 
ator mentioned above [14]. They tested this parallel-plate regener¬ 
ator using the same engine in which the stacked-screen 
regenerator had been tested. It has been reported that, using the 
parallel-plate regenerator and under similar test conditions, the 
engine delivered about 1750 W of acoustic power to its resonator 
[20], which was almost double the power delivered in the 
stacked-screen regenerator configuration [14]. This showed that 
the parallel-plate regenerator was clearly superior to the stacked 
screen one. The difference in regenerator performance was attrib¬ 
uted to a nonlinear resistance in the stacked-screen regenerator at 
high pressure amplitude, however, no further detailed analysis of 
this difference was presented. 

Yu et al. [21] carried out experimental investigations of the 
optimum transverse dimensions of the channel of the stacked- 
screen regenerator in a small-scale travelling-wave engine. The 
experimental results revealed that the ratio of 3 k lr h plays an impor¬ 
tant role in controlling the onset of the engine. Here the thermal 
penetration depth 3 k = y/2k/cop m c p (/< is the thermal conductivity, 
c p is the isobaric heat capacity per unit mass, p m is the mean gas 
density and co is the angular frequency). Using the lowest onset 
temperature difference as an indicator of regenerator performance, 
it has been shown experimentally that an optimum 3 k lr h was be¬ 
tween 4 and 7, when the hydraulic radius of the regenerator and 
the mean pressure vary in a large range (the investigations having 
been conducted for both helium and nitrogen). Furthermore, the 
experimental optimum 3 k lr h increases slightly as hydraulic radius 
decreases. In fact, decreasing the hydraulic radius increases the vis¬ 
cous resistance in the regenerator. However, the impact of these 
parameters of the regenerator on the optimum transverse dimen¬ 
sion was not investigated systematically. 

Some similar work on optimising regenerators has also been 
carried out for coolers. The non-dimensional parameter StPr 2l3 lf F 
was proposed to evaluate the regenerators with different geome¬ 
tries [22], where St is Stanton number, Pr is Prandtl number, and 
fp is a Fanning friction number. This parameter is a measure of 
the ratio of heat transfer to pressure loss. Mitchell et al. [23] pro¬ 
posed a “figure of merit” based on experimental data. It measures 
the heat transfer effectiveness per unit of flow resistance of the 
regenerators, and is expressed as follows, 



1 


M 


Pe , 4 Nk\ 
4 Nu ^ Pe ) 



where f D is the Darcy friction factor, Pe is Peclet number, Nu is the 
Nusselt number, and Nk is a measure of thermal dispersion. Based 
on this figure of merit, their experimental results showed that the 


etched stainless steel foil regenerator (similar to parallel-plates) 
was better than screens or fibers of 70% porosity, or packed spheres 
of 39% porosity. However, it should be noted that with a cooler, the 
thermal mass of the regenerators is much more important than that 
in the engines because of the higher heat capacity of the cold fluid 
and the lower heat capacity of the regenerator material at low tem¬ 
peratures. Furthermore, the local acoustic impedance of the regen¬ 
erator plays an important role in the travelling-wave engines. 
Therefore, it is necessary to include the effect of impedance for 
the comparison of regenerators. 

Generally, pressure and velocity amplitudes and their relative 
phase angles determine the energy fluxes through the regenera¬ 
tor. The impedance ultimately determines how large these fluxes 
are. The impedance values are generally determined by the 
acoustic network which consists of the acoustic compliance, 
inertance tube and the resistance. According to the literature 
survey above, it is clear that, unlike in the traditional Stirling en¬ 
gine, the role of the regenerator’s flow resistance in the travel¬ 
ling-wave thermoacoustic engine is complicated by the fact 
that the compact network needs to tailor both the impedance 
and the time phasing in the regenerator at the same time. It is 
still not very clear how the impedance affects the performance 
of the travelling-wave thermoacoustic engine in terms of the 
power output, onset temperature, and thermal efficiencies. 
Clearly, for different applications, different performance charac¬ 
teristics are preferred. Furthermore, the role of the flow resis¬ 
tance in the thermoacoustic regenerator also needs further 
study, together with its relationship with the regenerator 
impedance. 

This paper will focus on the regenerator of the travelling-wave 
engine, and will analyze the relationship between the acoustic 
impedance, flow resistance, viscous losses, heat relaxation loss 
and channel transverse dimensions. Much of the motivation be¬ 
hind the current work is driven by the interest as to whether 
there is a theoretical “optimum” of these parameters and/or the 
optimum of their combination for a specific application. Secondly, 
and more generally, the research attempts to develop a simple 
but effective methodology for the comparison of regenerators 
with different geometrical configurations or dimensions. In this 
work, the regenerators are studied in the travelling-wave for 
the ideal gas, based on the linear thermoacoustic theory. A pair 
of non-dimensional factors is obtained from the viewpoint of 
the “net” acoustic power production within the regenerator. 
These factors reflect the impact of the channel dimension, local 
acoustic impedance and the channel geometries on the power 
output capacity of the regenerators. Using this simple method, 
regenerators with different dimensions and channel geometries 
have been investigated numerically and compared with each 
other. 


2. Theoretical analysis 

2.1. Regular geometries of regenerators 

According to the linear thermoacoustic theory [24,25], the time- 
averaged acoustic power dE 2 produced in a length dx of the channel 
can be written in the complex notation in the general form as 


dE2 

dx 




dp, 

dx 


+ Pi 


dU{ 
~dx ’ 



where Ifl and p i are complex volumetric velocity and pressure, 
respectively. indicates a complex conjugate. Subscript 2 indi¬ 
cates that it is a second order quantity. Re[ ] denotes the real part 
of a complex number. Eq. (2) can be written as [25] 
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Viscous resistance per unit length of the channel, r v , thermal- 
relaxation conductance per unit length of the channel, l/r /o and 
the complex gain/attenuation constant for the volume flow rate, 
g, are defined as follows: 


r v 

1 

r k 


cop m Im[-/ t 


A 


g 


y - 1 coA g 


y P 


1 

Im[—/ fc ], 


m 


(4) 

(5) 


and 


g _ (fk fv) _1_ dT m 

(\-f v )(\-o)T m dx [ } 

f v and f k are given in [24] and [25] in detail, y, cr, p m , p m and T m are 
the ratio of specific heat capacities, Prandtl number, mean density, 
mean pressure and mean temperature of the working gas, respec¬ 
tively. A g is the cross-sectional area of gas channels in the regener¬ 
ator. Furthermore, A g = fA, where 0 and A are porosity and cross- 
sectional area of the regenerator, respectively. 

On the right hand side (RHS) of Eq. (3), the first two terms 
represent viscous and thermal-relaxation dissipation, respec¬ 
tively, which always consume acoustic power, regardless of the 
temperature gradient along the length of the regenerator. The 
third term denotes the acoustic power produced (or consumed) 
by the regenerator due to axial temperature gradient. It depends 
on the magnitude and direction of the axial temperature gradi¬ 
ent. In the engines, T m increases in the direction of positive 
acoustic power flow, so the third term denotes the acoustic 
power produced from heat energy. For the refrigerators, T m de¬ 
creases in the direction of positive acoustic power flow, in which 
case the third term means the acoustic power consumed to 
pump heat from the cold to the hot end of the regenerator. This 
paper focuses only on the analysis of regenerators in engines. 
Therefore, it will be more convenient to refer to dE 2 /dx as “net” 
time-averaged acoustic power production per unit length of the 
regenerator, because it is a net effect due to the acoustic power 
dissipation (the first two terms of the RHS of Eq. (3)) and the 
acoustic power production (the third term). 

For the travelling-wave, the phase angle between U 1 and p i, 
6(U lt p i) = 0. Hence 



Pm 

A 



Subscript “7W” refers to the travelling-wave, n is a real number 
and denotes a dimensionless local acoustic impedance factor. 
Substituting Eqs. (4)—(6) and (11) into Eq. (8), the following rela¬ 
tionship is obtained: 


^^2 1 | j j 1 2 J G)P m 

d^ = 2 m 


+ J-S=Re 





1 OjAg c t f Pm ^ 

n" m 


T m dx L(1 -/»)(! - cr)J V A? 



A 


g 


( 12 ) 


For ideal gases, the speed of sound a = ^/yM\/ is the gas 
constant per unit mass), and the mean pressure p m = p m MT m . In 
addition, a = coy/2n, where 2 is the wavelength. Thus, Eq. (12) 
can be simplified to: 


dEj 

dx 



-71 


Im [~fv 
|1 -fv\ 



(fk -fv) 

dT m /dx \ 

[(1 -/„)(!- <7)J 

T m /X J 


nn 2 (y- l)Im[—/ fe ] 


.(13) 


The factor in the curly Brackets is dimensionless and in the cur- 


It can be found that H 


|l/i 1 2 has the same dimension as that of 


dx 


rent analysis is expressed as follows: 


'F 


U 


-71 


|1 -fv 


- 7trr(y - l)Im[-/ k ] 


(fk -fv) 


+ -Re 

+ 2 [(I-/,)(!-<r)J 


0 


VT: 


(14) 


where 6> V r is a dimensionless temperature gradient defined as 



dT m /dx 

~^m/T 



Therefore, dE/dx is proportional to Wu if ^ gi ven - 

On the other hand, substituting Eq. (11) into Eq. (13) yields 


dE 2 

dx 



V mh/l 

n 2 


|1 -fv\ 


7t(y- l)Im[—/ k ] 


(fk ~fv) 

dT m /dx\ 

[(1 -/„)(! - <T)J 

T m /X j 


(16) 


2 

Similarly, g ^lso has the same dimension as that of The 
factor in the curly brackets is dimensionless and can be expressed 
as follows: 


Re[p 1 l/ 1 ] = |p 1 ||l/ 1 |, Im[p 1 H 1 ]=0. 

Substituting Eq. (7) into Eq. (3), the following is obtained: 

_ 1 ... af 1 IP! | 2 p w li 

In Eq. (8) 

|Pil/|Ui| = |Pi/U,|, 


(7) 

V P = \n 11 

y n 2 |i 

(8) 

+ } Re 
2 n 


( 9 ) 


1 ~fv\ 


n(y- l)Im[—/ fc ] 


(fk -fv) 


L(i -f»)0 ~ ff )J 


e 
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(17) 


Here, we call Wu and Wp “power factors” because they reflect the 
power output capacity of the regenerators. One can find that 


W 


V 


u 


p 




(18) 


where pi/lfi is the local acoustic impedance Z of the acoustic wave. 
For a lossless planar travelling-wave, Z is defined as [24] 


y _ Pi _ Pm^ 

~ur A g 



where a is the speed of sound. However, for the practical travel¬ 
ling-wave engines, the local acoustic impedance should be much 
bigger than p m a!A g to decrease the large viscous losses resulting 
from high acoustic velocities [9,12]. Here, for convenience it is ex¬ 
pressed as 


According to the derivation above and Ref. [24], it can be found 
that (or Wy) is a function of f k ,f v , n, cr, and y. For regular geom¬ 
etries, f v and f k are analytical functions of S k lr h • cr and y are given 
when the working gas is given. Assuming the temperature distri¬ 
bution is linear along the length of the regenerator, the mean value 
of 6>vr can be estimated from Eq. (15) by using the total length of 
the regenerator, the temperature difference between the two ends 
and the mean temperature of the gas in the middle of the stack 
(T m ). Therefore, the impact of Sjjr h and n on dE/dx can be studied 
numerically. 
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Compared with Eq. ( 8 ), it can be found that, in the RHS of Eqs. (14) 
and (17), the first and second terms denote dissipation by viscosity 
and thermal relaxation, respectively. The third term represents the 
acoustic power produced from heat energy. Qualitatively it can be 
concluded that, the acoustic impedance reflects the proportion 
among the viscous dissipation, thermal-relaxation dissipation, and 
the new acoustic power production from heat energy. According to 
Eqs. (14) and (17), the weighting of the third term is much bigger 
than that of the first term, because <9 V r > 1 and n > 1. Therefore, 
one can see that increasing volumetric velocity increases the first 
term, which means viscous dissipation, however, it also increases 
the third term which corresponds to the acoustic power production. 

However, it should be noted that, in practical travelling-wave 
engines, |Lfi|, |pi| and n are determined by the network of all the 
acoustic components which is often quite complex. One cannot 
arbitrarily increase or decrease n whilst keeping |Lfi| or |pi| con¬ 
stant. Nevertheless, one could focus on the analysis of how the 
transverse dimension and impedance impact the acoustic power 
flow in the regenerators. Therefore, Eqs. (14) and (17) will give a 
useful insight into the regenerator’s power output capacity. 

2.2. Stacked-screen regenerators 

For most practical travelling-wave thermoacoustic systems, 
the regenerator comprises of layers of stacked screens (wire- 
meshes). It is possible to perform the analysis similar to that 
outlined in Section 2.1, providing appropriate analytical expres¬ 
sions can be found. Swift and Ward [27] proposed a set of ther¬ 
moacoustic equations for stacked-screen regenerator with the 
method of simple harmonic analysis. Based on their equations, 
the analysis of the acoustic power dissipation and production 
in the stacked-screen regenerator can be developed further. 
Substituting their continuity and momentum equation [27] into 
Eq. (3) yields 




ncj f3 ,<\ /ci(0) c 2 ((p)Re i 

n~ + ^r 


nn 2 (y - 1 ) 


x Im 
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Here, subscript “SCR” indicates a stacked-screen regenerator. An 
analogous equation to (17) has the following form: 
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Here, 'Fu,scr and W PtS cr are the pair of “power factors” for stacked- 
screen regenerators. Similarly, it can be found that 

'Pp,scr = 'Puscr/h 2 - ( 25 ) 

According to the derivation above and Ref. [27], it can be found 
that Wp'SCR (or Y UrSCR ) is a function of 3 k lr h , n, </>, Re 1 , 0 VT , cr, 7 , e s 
and 0 T . Furthermore, for practical travelling thermoacoustic engine, 
the impact of e s and 0 T on ^u.scr can be neglected [27]. cr and y are 
given when the working gas is given. Therefore, the impact of 3jJr h 
and n on dE 2 /dx can be studied when </>, Re t and <9 V r are given. 


3. Numerical analysis 
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In Eq. (19), ( ) denotes a local spatial average. p m , y m and \x m 
are mean thermal expansion coefficient, mean ratio of specific 
heats, and mean dynamic viscosity of the gas, respectively. b(4> 
), Ci( 0 ), c 2 (0), gc(Rei), and g v (Rei) are factors resulting from the 
fitting of data of Kays and London [27,28]. 6 P and 0 T are defined 
as 0 p = 0(Pi,<Ui» and 0 T = 0«Ui), <T) Ui ,), respectively, e s = 4’PmC p ,ml 
(1 - (t>)ps,mC s ,m (subscript s means solid) gives the ratio of gas 
heat capacity to solid heat capacity in the stacked-screen regen¬ 
erator. s h is defined as s h = 8ir^/h(0)cr 1 / 3 ^. Reynolds number 
Rei =4|(Ui)|rhp m //i m is a positive real number. More detailed 
description of these variables and definitions is given in Ref. 
[27]. 

By analogy with Eq. (11), the effective local specific acoustic 
impedance can be introduced as follows: 


z e =p 1 /(uj) =np m a. (20) 

Then, assuming a pure travelling-wave and ideal gas, one can 
obtain Eqs. (21)-(23) for stacked-screen regenerators through sim¬ 
ilar derivations as those for the regular geometries 


dE 2 

dx 


( 27 T^ Pm K Ul >| 2 ) lFu ' SCR ^ 



Theoretically, both W p and could be used for the analysis of 
the regenerators. However, in practice, it is more convenient to 
measure \p^\ than |Lfi|. Therefore, W p and W p; S cr are used in the 
analysis below, where numerical calculations for parallel-plate 
and stacked-screen regenerators have been performed. The dimen¬ 
sionless transverse dimension of the channel, 3 k lr h , was used for 
the analysis. All calculations were performed for helium as the 
working fluid (cr = 2/3, y = 5/3). 

According to Eqs. (17) and (24), <9 V r is a required value for the 
numerical analysis. As shown in Eq. (15), <9 V r can be estimated by 
using the length of the regenerator instead of dx, and the temper¬ 
ature difference between the two ends instead of dT m with an addi¬ 
tional assumption of a linear temperature gradient along the 
regenerator. The wavelength is decided by the operating frequency 
(usually several resonator lengths in practical systems). Therefore, 
one can obtain the approximate mean value of 0 VT in the regener¬ 
ator. For example in Backhaus and Swift’s travelling-wave engine 
[7], at its point of the highest power production, the hot and cold 
end temperatures of the regenerator are 725 °C and 40 °C, respec¬ 
tively. The length of the regenerator is 7.3 cm. The operating fre¬ 
quency is 84 Hz with helium at a mean pressure of 31 bars as 
the working gas. Based on this information, mean value of <9 V r is 
estimated as about 172 in the regenerator. Here, for the conve¬ 
nience of further comparison with these experimental results, 
©vt= 172 is set for the calculations. 

3.1. Parallel-plate regenerator 

Fig. 2 shows the calculated Wp for the parallel-plate regenerator 
when 3 k lr h varies, for several values of n treated as a parameter. 
There is a maximum value of *R P for each n when 3 k lr h varies, 
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Fig. 2. Parallel-plate regenerator in the travelling-wave: the calculated value of W p 
versus 5 k /r h for different local acoustic impedances. For each n, there is a maximum 
value of W p when d k /r h varies. 



n 

Fig. 3. The relationship of ( W p ) max versus n for both the parallel-plate and stacked- 
screen regenerator. The two plotted points schematically show the operating states 
of the two regenerators tested by Backhaus and Swift [7,20]. 


designated as ( W p ) max . This corresponds to the optimum value of S k l 
r h , denoted as ( S k lr h ) opti . On the other hand, one can find that the opti¬ 
mal range of S k lr h becomes broader as n increases. For example for 
n = 25 and n = 35, *F P increases sharply when S k /r h < 2.5, but for S k l 
r h > 3 the dependence is rather weak. This means that in practical 
travelling-wave thermoacoustic engines, a broad range of S k lr h 
may be suitable for regenerators. Furthermore, *F P decreases rapidly 
as n increases if S k lr h is fixed. Because !F P is related to the net time- 
averaged acoustic power production of the unit length of the regen¬ 
erator, this figure reveals that both the local acoustic impedance and 
channel dimension significantly influence the acoustic power out¬ 
put of regenerators in the travelling-wave engines. 

In addition, Fig. 2 can be interpreted as a graph of (W p ) max and 
( Sklrh) op ti as functions of n (see solid lines in Figs. 3 and 4). In Fig. 3, 
the solid line shows that ( l F p ) max is nearly in inverse proportion to 
n. This agrees with the analysis in Section 2.1 : the third term in Eq. 
(17), which has the biggest weighting, is indeed inversely propor¬ 
tional to n. Flence, if g is given, the acoustic power output of 
the regenerator is in practice almost inversely proportional to 
impedance. 

Fig. 4 shows the relationship of ( $ k lr h ) opti versus n. It can be seen 
that, ( S k lr h ) opti slightly increases as n increases. Clearly, the opti¬ 
mum dimensionless channel dimension in the regenerator exists; 
furthermore, it depends on the local impedance. As the impedance 
increases, the corresponding optimum ( $ k lr h ) opti increases, i.e. the 
optimum traverse dimension of channel decreases. Theoretically, 
it is possible to obtain an optimum combination between the 
transverse dimension r h of channel and impedance, which would 
lead to the highest acoustic power output. 

3.2. Stacked-screen regenerator 

Eq. (24) is used for the analysis of the stacked-screen regenera¬ 
tors. Compared with those for parallel-plate regenerator, two 
additional parameters, 0 and Re A should be provided for the calcu¬ 
lations. To compare with the results of Backhaus and Swift’s 
travelling-wave engine (mean pressure of 30 bars) with the 
stacked-screen regenerator [7], <9 V r, 4> and Re\ are set according 
to its highest power operating point. In this case, <9 V r is also esti¬ 
mated at about 172, and the measured porosity 0 = 0.72. One can 
estimate the mean value of Re i at about 50 in the regenerator 
based on DeltaE [29] simulations. 


6 



CO 


w 2-'" 

1 _ -Parallel-plate, © V1 =172 

.Stacked-screen, 0 VT =172, /te =50. <|> =0.72 
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Fig. 4. The calculated relationship between {SiJr h ) opti and n for the tested parallel- 
plate and stacked-screen regenerators. 



Fig. 5. Stacked-screen regenerator in the travelling-wave: the calculated values of 
*p. scr versus S k lr h for different local acoustic impedances. For each n, there is a 
maximum value of W PiS cr when S k /r h varies. 
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The calculated results are quite similar to those for the parallel- 
plate regenerator and shown in Figs. 3-5. Firstly, Fig. 5 shows the 
calculated Wp^cR for the selected stacked-screen regenerator when 
S k /r h varies, for several values of n. The optimal points in Fig. 5 are 
also shown as dotted lines in Figs. 3 and 4, respectively. 

From Figs. 3-5, it can be found that the local acoustic imped¬ 
ance and channel dimensions influence the power output of 
stacked-screen regenerator almost in the same way as those for 
the parallel-plate regenerator. The method developed for the regu¬ 
lar geometrical configuration is applicable to the stacked-screen 
regenerators. 

3.3. Comparisons with experimental results 

As already mentioned, the regenerators discussed above have 
the same dimensions as those studied by Backhaus and Swift 

[7.20] . Their experimental results indicate that, under nearly the 
same operating conditions, the maximum acoustic power output 
of the travelling-wave engine with the parallel-plate regenerator 
is almost double of that for the engine with the stacked-screen 
regenerator. Flowever, according to Fig. 3, one can see that the dif¬ 
ference between (*F p ) m ax and (^.sck) max is very small (about 4%) 
for each given n. The reason behind such a large difference in 
power output needs to be identified. 

For further analysis, DeltaE [29] code is used to simulate the 
whole engine of the two cases tested by Backhaus and Swift 

[7.20] , based on the published experimental results. Dividing the 
regenerator into 16 segments, the distribution of pressure ampli¬ 
tudes |pi| and volumetric velocity |Lfi| are obtained and shown in 
Fig. 6a and b, respectively. In Fig. 6a, it can be found that |pi| de¬ 
creases along the length of the regenerators from the cold to hot 
end. Furthermore, the pressure amplitude drops |Api| are about 
0.3 and 0.15 bar in the tested stacked-screen and parallel-plate 
regenerators, respectively. In Fig. 6b, it can be seen that the volu¬ 
metric velocity |Lh| increases along the length of the regenerators 
from the cold to hot end. Flere, one can find that the temperature 
gradient has a large effect on the volumetric velocity. Furthermore, 

| (Ji| through the parallel-plate regenerator is about 1.5 times larger 
than that through the tested stacked-screen regenerator at each 
cross section. Based on the simulated results shown in Fig. 6a 
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and b, one can calculate the approximate dimensionless imped¬ 
ance factor n in the two tested regenerators using Eqs. (11) and 
(20), as shown in Fig. 6c. It can be seen that n decreases along 
the length of the regenerator from the cold to the hot end. Further¬ 
more, the mean value of dimensionless impedance factor n in the 
stacked-screen regenerator is about 1.7 times of that in the tested 
parallel regenerator. 

Using the above simulation results, the viscous dissipation in 
the tested stacked-screen regenerator is estimated as about 40 W 
more than that in the parallel-plate regenerator. Therefore, the vis¬ 
cous dissipation difference could be partly a reason behind the 
large output difference. However, this is still not large enough to 
explain the experimental differences. 

Fig. 6d shows the obtained phase angle between Ui and pi 
through the regenerators. It can be seen that, in the two regenera¬ 
tors, the phase angles 0(lfi,pi) are a bit different near the cold end, 
but quite similar near the hot end of the regenerators. This differ¬ 
ence is mainly due to the different resistance in the two regenera¬ 
tors. Given that the majority of the acoustic power is produced 
near the hot end of the regenerators, the difference of the phase an¬ 
gle shown in Fig. 6d will also account for some of the difference in 
the acoustic power outputs of these test cases, but clearly it is still 
not the key reason. 

According to Fig. 6c, the estimated values of n close to the mid¬ 
dle of the regenerators are about 14 and 25 in the parallel-plate 
and stacked-screen regenerators, respectively. Going back to 
Fig. 3, the two tested regenerators by Backhaus and Swift [7,20] 
can be schematically shown by the two plotted points. One can 
find the corresponding (Wp ) max ~ 5.9 and (‘f'p.scfOmax ~ 3.1, which 
indicates that the possible maximum capacity of acoustic power 
output for the parallel-plate regenerator is about twice of that from 
the stacked-screen regenerator. It is quite interesting that the ob¬ 
tained result agrees with the experimental results qualitatively 

[7,20], although this kind of analysis is rather approximate, be¬ 
cause it is based on an ideal gas and a pure travelling-wave. There¬ 
fore, it can be qualitatively concluded that the difference of the 
impedances is the possible key reason behind the large difference 
in power outputs for the two tested cases. 

3.4. The role of the viscous resistance in the regenerator 


At this point, the relationship between flow resistance and 
impedance can be discussed further. In the same notation as Eq. 
(4), the viscous resistance per unit length of stacked-screen regen¬ 
erator is [25,27] 


f v,SCR 


Vm ( c iW c 2 (<ft)ReA 
r 2 h 4>A\ 8 3n )' 


(26) 


Substituting the viscous penetration depth S„ = x /2fi ln /oip m 
into Eq. (26) to cancel /u m leads to 


f v,SCR 


1 cop m fd v \ 2 fc^) c 2 (<p)ReA 

2 ,j>A \r h J V 8 3n )' 


(27) 


Based on the parameters used in Sections 3.1 and 3.2, one can 
estimate the mean resistance along the length of the two regener¬ 
ators by using Eqs. (4) and (27). The mean value of r v>sc Wrv in the 
regenerators is about 2.8. So the viscous resistance of the 
stacked-screen regenerator is substantially bigger than that of par¬ 
allel-plate regenerator. 

On the other hand, neglecting the compliance effect due to 
regenerator’s porosity, the pressure amplitude drop through the 
regenerator [7] 


Fig. 6. Results of DeltaE simulations, based on the experimental data of Backhaus 
and Swift [7,20]. Distributions of |pt|, |lh|, n and 6(11-1, Pi) through the two tested 
regenerators. 
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Here, t = T h IT c , the subscripts h and c mean the hot and cold end 
of the regenerator, and R v is the viscous resistance of the regen¬ 
erator. Therefore, according to Eq. (28), it can be seen that higher 
viscous resistance R v (or r v ) leads to a smaller and thus will 
lead to a high impedance when the pressure amplitude drop Ap i 
is given. This analysis is consistent with the simulation results in 
Fig. 6 . 

As discussed in the introduction, the thermoacoustic core in 
such travelling-wave engines is an acoustic power amplifier. This 
amplification of acoustic power can be approximately expressed 
as, 

E„ut = -cEin, (29) 


where E in and E out are the acoustic power that runs into the ambient 
end of the regenerator, and the acoustic power runs out of the hot 
end of the regenerator, respectively. Therefore, the total net acous¬ 
tic power produced by the regenerator can be defined as [7] 

E ne t = (T-l)Ein. (30) 


For a given pressure amplitude at the ambient end of the regen¬ 
erator pi >in , W in can be approximately expressed as, 



|Pl,inl 

2Re[Zj„] ’ 



Assuming a travelling-wave phase at the ambient end of the 
regenerator, substituting Eq. (11) into Eq. (31) leads to (cf. [19]) 


4IPl,.n! 2 

2 np m a ' 



By combining Eqs. (32) and (30), the relationship between the 
impedance and the net acoustic power production can be obtained 
as 



(T - l)4|p 1|in 
2 np m a 



Eq. (33) indicates that, if the operating conditions and the pres¬ 
sure amplitude are given, the net power production in the regener¬ 
ator is almost inversely proportional to the impedance factor n, 
which agrees with Eqs. (17) and (24), and the results shown in 
Fig. 3. 

Qualitatively, it can be concluded that a higher viscous resis¬ 
tance of the stacked-screen regenerator leads to a smaller volumet¬ 
ric velocity, and consequently, also a higher impedance, in the 
tested stacked-screen regenerator [7] compared with the paral¬ 
lel-plate regenerator [20]. A smaller volumetric velocity and a 
higher viscous dissipation both result in a smaller acoustic power 
output for the travelling-wave engine with the stacked-screen 
regenerator [7]. But, the smaller volumetric velocity (i.e. higher 
impedance as shown in Fig. 3) is the key reason because it de¬ 
creases the power input fed back to the cold end of the stacked- 
screen regenerator [7] compared with the parallel-plate one [20]. 

It should be noted that the current analysis does not account for 
any nonlinear viscous resistance in the stacked-screen regenerator 
at high acoustic amplitudes. If this nonlinear resistance were con¬ 
sidered, r v<S cR^v would increase. Thus, the volumetric velocity 
would decrease in the stacked-screen regenerator, and subse¬ 
quently the impedance would increase. As a result, the predicted 
difference of power output capacity as shown in Fig. 3 would be 
higher. 


4. Discussion and conclusions 

Based on the assumptions of the pure travelling-wave and ideal 
gas, a pair of dimensionless factors, Wu and W p , have been obtained 
to evaluate the acoustic power output capacity of regenerators in 


the travelling-wave thermoacoustic engines. These factors can 
clearly reveal the relationship between the flow resistance, imped¬ 
ance, and channel dimensions with the acoustic power output 
capacity. 

It is shown that the local acoustic impedance of the regenera¬ 
tors essentially represents the proportion between the acoustic 
power produced from the heat energy through the thermoacoustic 
processes and the acoustic power dissipated by viscous and ther¬ 
mal-relaxation effects in the regenerators. The net acoustic power 
output of the regenerator is determined by its impedance. This is 
because the impedance determines the amount of power fed into 
the ambient end of the regenerator. For a given pressure ampli¬ 
tude, a higher impedance leads to a lower acoustic power input. 
Therefore, the local acoustic impedance has to be taken into ac¬ 
count for any evaluation of regenerators, as well as the comparison 
between different geometrical configurations or dimensions. 

For a given pressure amplitude, |pi|, increasing the impedance 
will decrease the volumetric velocity, and thus will decrease the 
acoustic power output of the engine. Therefore, relatively low 
impedance is preferred for the travelling-wave thermoacoustic en¬ 
gines to obtain a high acoustic power output. Certainly, higher vol¬ 
umetric velocity will cause higher dissipation if the resistance is 
not reduced, and then degrade the thermodynamic efficiency. 
There should be a compromise between high efficiency and high 
power output according to the requirements of different 
applications. 

The resistance of the regenerator plays the key role in determin¬ 
ing the regenerator’s impedance. For a given pressure amplitude, 
the higher resistance will decrease the volumetric velocity through 
the regenerator, which consequently causes the high impedance 
of the regenerator. This is the key reason of the large difference 
of power output for the tested stacked screen [7] and parallel-plate 
[20] regenerators. In any case, decreasing the flow resistance is 
most crucial in the practical design of regenerators for travelling- 
wave thermoacoustic engines. Practically, one can reduce the 
regenerator’s resistance by using regular configurations such as 
parallel-plate to achieve relatively low resistance, reducing the 
length of the regenerator, or increasing the cross-sectional area. 

Using the same methodology as in Section 3.1, one can perform 
a similar numerical analysis for pin-array and circular-pores regen¬ 
erators. Obtaining similar results to those shown in Figs. 2-4 is rel¬ 
atively straightforward for these two geometrical configurations 
and therefore has been omitted in this paper. Compared with par¬ 
allel-plate regenerator with the same porosity, it is found that for 
given impedances these regenerators have almost the same possi¬ 
ble maximum acoustic power output capacity as those of parallel- 
plate regenerator. In other words, from the point of view of Rott’s 
function f k , the geometrical configuration is not a sensitive factor to 
regenerator’s performance for a given impedance. This is mainly 
due to the fact that for all these geometrical configurations Re\f k ] 
approaches 1 when d k lr h is high [25]. However, the regenerator’s 
resistance changes this picture. As shown in Fig. 3, for given 
impedance, the difference in W p between the stacked screen and 
the parallel-plate regenerator is quite small. But different imped¬ 
ance due to different resistance caused large differences in power 
output. Therefore, it can be inferred that, regardless of the geomet¬ 
ric configuration, the regenerator with lower resistance will have a 
better performance when the other parameters (i.e. porosity, chan¬ 
nel dimensions, material, etc.) are the same. 

This mechanism is different from that behind the large effi¬ 
ciency difference due to different geometries of stacks in stand¬ 
ing-wave engines. It was reported that, among the three regular 
geometric configurations mentioned above, the pin-array stack is 
the most efficient one, parallel-plate is the moderate, and the cir¬ 
cular-pore stack is the worst one [26,30]. For the stacks, the perfor¬ 
mance difference is mainly due to the fact that \m[f k ] is quite 
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different when 3 k lr h ~ 1. In contrast to the regenerators, the stacks’ 
impedance is determined by its location in the standing-wave res¬ 
onator rather than its resistance. 
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